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ON  THE  STRl’lTUKE  OK  LAMINAE  DIFFUSION  FLAMES 


AMAHLK  LIN  A  N 


A  it  S  T  It  A  C  T 

Thu  structure  of  1  a n i i n;ir  diffusion  flumes  is  analyzed  in  tin?  limiting  rasa  of  large,  although  finite,  reaction 
rates. 

it  is  shown  that  the  ehemiea!  reaction  takes  place  only  in  a  very  thin  region  or  “Chemical  boundary  layer* 
where  convection  effects  may  he  nejrleeted.  Then  the  temperature  and  mass  fraction  distributions  within  the 
reaction  zone  are  obtained  analy.ically, 

Tlie  flame  position,  rates  of  fuel  consumption,  and  temperature  and  concentration  distributions  outside  of 
the  reaction  zone  may  he  obtained  by  using  the  assumption  of  infinite  reaction  rates. 

Kor  large  Reynolds  numbers  mixing  and  combustion  take  place  in  boundary  layers  aud  free  mixing  layers. 
And  again  analytical  solutions  are  obtained  for  the  temperature  and  mass  fraction  distributions  outside  of  the 
reaction  zone. 


X  <>  M  H  NCI,  A  T  II  R  K 

Tltt'  following  is  tt  list  of  (lie  most  important  symbols  used  in  this  paper. 

A„  I’aramotor  given  by  [44 1.  that  measures  the  deviations  from  the  Ilurko-Sehuniann 
solution. 

r,,  Speeifie  heat  tit  eonsttmt  prt'ssure. 

It  Diffusion  coefficient. 

K  Activation  energy  of  the  ehemiettl  reaction. 

/  Dimensionless  stream  function. 

K,  Mttss  fraction  of  species  /. 

L  Some  overall  eharaeteristie  length. 

M  Mean  molecular  mass. 

m  Mttss  rate  of  fuel  eomsumption  per  unit  flame  surface, 

l’r  ITandtl  number. 

/i  l’ivssiire. 

(/  ileal  released  per  unit  mttss  of  fuel. 

I!  Universal  gas  constant. 

Re  Reynolds  number. 

Se  Schmidt  number. 

T  Temperature. 

T,  Adiabatic  flame  temperature  given  by  |24). 

Tr  Tempertiture  tit  the  ideal  flattie  surfttee. 

/,  Characteristic  chemical  time  defined  by  [15], 
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tm  Characteristic  mixing  time,  4=  Si/ IV 

U  Characteristic  overall  velocity. 

it,  v  Velocity  components  in  boundary  layer  coordinates. 

v  Velocity  vector. 

v,ij  Diffusion  velocity  of  species  i. 

H'j  Mass  production  rale,  per  unit  volume,  of  species  L 
X,  V  Mixing  boundary  layer  coordinates, 

a'  Position  vector. 

x,  n  Chemical  boundary  layer  coordinates. 

•I 

!h  =  ](?/  !>..)<///• 

0 

=  .'/i/'i- 

Universal  function  giving  the  temperature  distribution  within  the  reaction  zone.  De¬ 
fined  by  |:l!>). 

Mixing  length.  o„  D„/mi. 

\  Characteristic  thickness  of  (lie  reaction  zone  given  by  [41], 
r,  Dimensionless  distance  normal  to  the  mixing  layer. 

»  Non-dimensional  temperature.  9  =  (T  —  T0)/(T. — '!'„)• 

o„  =  k/khv  t„). 

9,1  —  1  rt/(  1  /  l  0  )- 

K  =  (Tr— T„)/(Tr  T„). 

). 

H  Viscosity  coefficient, 

v  Stoichiometric  ratio  species  2  to  1. 

;  Non-dimensional  distance  along  the  mixing  layer. 

|>  I  tensity. 

3  Stress  tensor. 

SrnscmiTs: 

1,2,;!  Indicate  fuel,  oxidizer  and  products  respectively. 

Indicates  conditions  at  the  fuel  exit. 

„  Indicates  conditions  on  the  oxidizer  side  of  the  flame,  far  front  tin'  flame. 

/'  Indicates  conditions  at  the  flame  surface  for  infinite  reaction  rates. 

The  asterisk  is  used  for  the  non-dimensional  variables  introduced  in  section  II.  (i. 


1  .  In  t  n  o  it  r  ct  t  <>  x 

I  til  fusion  flames  are  obtained  when  the  reacting  species  are  initially  separated.  ( 'out bas¬ 
tion  and  mixing  takes  place  simultaneously. 

In  these  flames  the  reaction  zone  separates  the  two  reacting  species  which  diffuse,  through 
inert  gases  and  combustion  products,  from  each  side  towards  the  flame. 

Tiie  reacting  species  burn  very  rapidly  as  they  reach  the  reaction  zone;  thereby  the 
combustion  velocity  is  generally  conditioned  to  the  aecesihility  of  the  species  to  the  reaction 


zono:  or  in  othor  words,  to  their  facility  to  diffuse  across  the  inert  gases  and  combustion 
products. 

It  seems  that  we  can  arrive  at  a  description  of  some  of  the  most  important  features  of 
diffusion  flames  by  using  the  assumption,  first  introduced  by  Burke  and  Schumann  (1),  of 
infinitely  fast  reaction  rates.  Then  the  actual  zones  of  combustion  become  infinitely  thin,  and 
the  mixing  process  alone  becomes  responsible  for  the  rate  of  burning  and  for  flame  location 
and  size. 

Iturke  and  Schumann  have  successfully  used  their  assumption  for  the  calculation  of  the 
shape  and  length  of  the  laminar  diffusion  flame  formed  when  a  fuel  jet  discharges  within  a 
tube.  In  this  lube  an  air  stream  moves  with  the  same  velocity  as  the  fuel  jet.  The  same  as¬ 
sumption  has  been  utilized  by  llottel  and  Hawthorne  (2),  Wold,  (iazley  and  Kapp  (:!),  \>gi 
and  Saji  (4).  and  Harr  (ii),  for  the  prediction  of  the  length  of  open  flames,  both  laminar  and 
turbulent.  Through  rudimentary  approximations  they  obtain  an  expression  for  the  flamo 
length  containing  an  unknown  function:  this  they  determine  empirically  from  the  results  of 
their  experiments.  Fay  (<»)  has  calculated,  by  using  Burke  Schumann  assumption,  the  shape 
and  characteristics  of  the  laminar  diffusion  flame  obtained  when  a  fuel  jet  discharges  into 
the  open  atmosphere. 

The  infinite  reaction  rale  assumption  has  also  been  utilized  (7).  (8)  for  the  study  of  diffu¬ 
sion  flames  in  boundary  layers. 

In  addition,  an  extensive  literature  exists  on  the  application  of  the  assumption  to  fuel 
droplet  combustion. 

The  Iturkc  Schumann  assumption  eliminates  chemical  kinetics  from  the  process,  simpli¬ 
fying  the  governing  cipialions  and  their  solution.  I  Iowever,  this  solution  does  not  provide  the 
criterion  for  the  extinction  of  the  flame,  or  for  the  validity  of  the  assumption  and  solution. 

Zeldovich  (it)  has  taken  into  consideration  the  finite  thickness  of  the  reaction  zone  to 
explain  the  blowing-off  phenomenon.  Similar  studies  have  been  performed  by  Spalding  (10), 
(I  I).  (12)  with  the  purpose  of  relating  the  fuel  consumption  rate  per  unit  area  at  extinction 
and  the  fuel  consumption  rate  per  unit  area  in  a  preniixcd  flame. 

For  a  general  description  of  the  diffusion  flames  see,  for  example,  the  review  papers  by 
Barr  (It))  and  Wold  and  Shipment!  (14).  where  data  and  bibliography  on  the  subject  can  be 
found. 

We  aim  in  this  work  to  show  the  effects  of  finite  chemical  reaction  rates  on  the  structure 
of  laminar  diffusion  flames.  In  order  to  do  so,  we  will  study  certain  limiting  eases,  in  which 
simple  analytical  solutions  can  be  obtained.  We  will  limit  ourselves  to  the  study  of  one  step 
chemical  reactions  in  which  the  forward  reaction  is  dominant. 

We  shall  show  that  for  large  reaction  rates  the  chemical  reaction  takes  place  only  in  a 
very  thin  region  or  .chemical  boundary  layer-.  This  has  already  been  shown  (15)  in  the  simple 
case  of  the  mixing  and  combustion  of  two  parallel  streams  of  fuel  and  oxidizer  moving  with 
the  same  velocity.  There  convection  effects  may  lie  neglected  compared  with  the  much  more 
important  diffusion  conduction  and  chemical  reaction  effects.  The  governing  e<|imtions 
reduce  in  this  cast'  to  ordinary  differential  equations.  The  kinetics  of  the  reaction  appears 
in  the  solution:  but  the  temperatures  arc  close  to  ilic  adiabatic  flame  temperature,  and  in  this 
range  of  temperatures  the  concept  of  an  overall  kinetic  scheme  has  been  found  by  Levy  and 
Weinberg  ( Hi)  to  be  valid. 

The  solution  with  the  assumption  of  infinite  reaction  rates  (which  we  shall  cal!  the  Burke- 
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Schumann  solution)  represents  the  true  solution  outside  of  the  reaction  zone.  It  may  also  be 
used  to  calculate  the  flame  position  and  fuel  consumption  per  unit  flame  area. 

If  the  Reynolds  number,  based  on  some  overall  characteristic  length,  is  large,  mixing  and 
combustion  will  take  place  only  in  a  very  thin  region  or  mixing  layer,  where  boundary  layer 
approximations  may  lie  used  (17),  (18),  (1!)). 


INVISCID  REGION 


MIXING _ /  |  OXIDIZER  SIDE 

REGION 


I'i”.  1. — I'iffuM.iii  KIiiiim*  Strin'tur*- 


!>  - 


The  mixing  layer  location  and  general  flow  characteristics  outside,  of  the  mixing  layer 
may  he  determined  by  using  the  inviseid  flow  equations.  However,  we  must  allow  for  the 
existence  of  discontinuities  in  the  velocity,  density,  temperature,  and  mass  fraction  distribu¬ 
tions  within  the  flow  field. 

In  Figure  1  the  temperature  and  mass  fraction  distributions,  as  obtained  by  different 
limiting  assumptions,  are  schematically  represented. 

II.  (iKN'KKAI.  KOt’ATIONS 

We  shall  begin  by  writing  the  general  ('((nations  governing  the  steady  laminar  flow  of  a 
reacting  gas  mixture  (20),  (21 ),  (22).  We  will  use  the  assumption  that  the  fluid  may  he  consider¬ 
ed  as  a  continuous  medium  formed  by  a  mixture  of  perfect  gases. 

Only  three  species  will  he  considered:  Fuel,  oxidizer,  and  products.  For  the  sake  of 
simplicity,  any  inert  species  present  will  be  considered  as  products. 

besides  the  usual  dependent  variables  of  ordinary  fluid  mechanics,  i.c.  velocity  v,  pres¬ 
sure  ;j,  density  p,  and  temperature  T,  three  new  variables,  the  mass  fractions  of  the  reactant 
species,  enter.  Therefore,  three  new  equations,  stating  the  mass  conservation  law  for  each  of 
the  species,  must  he  added  to  the  fundamental  equations  of  fluid  mechanics.  In  addition,  the 
relations  between  the  transport  parameters  and  mass  fractions,  temperature,  and  pressure  of 
the  mixture  will  he  required. 

We  shall  use  subscript  1  for  fuel.  2  for  oxidizer,  and  H  for  the  products.  The  mass  frac¬ 
tions  of  species  i  will  he  written 

K,  ii, iii. 

The  three  mass  fractions  obviously  satisfy  the  relation 

K,  |  K,  i  K:,  =  1.  [1] 


1 1.1.  Equation  of  Stair. 

If  the  fluid  is  considered  as  a  mixture  of  perfect  gases  the  equation  of  stale  is  as  follows 

l>  u  ( I  (M )  T  [2] 

where  It  is  the  universal  constant  of  the  gases,  and 


M  Vk  M. 


m 


is  the  mean  molecular  mass. 

lit  nrilri*  Ui  simplify  (hr  \v<*  will  um* 

proximal  ion  is  justified  when  the  molecular  masses 
when  the  reactants  are  very  dilute.  Then  M  M  ,. 


a  mean  constant  value  for  M.  This  au- 
of  t|ii»  him*  iiot  Very  <!i jf<»r<*n{  or 


In  any  ease  the  results  will  not  he  essentially  changed  by  considering  M  as  variable. 


II.2.  Equation  of  Continuiti)  for  the  Mixture. 


This  simply  states  the  law  of  mass  conservation 

V(pe)=0.  [4] 

II. U.  Equations  of  Mass  Conserration  for  the.  Species. 

These  state  that  the  mass  tpiuutity  of  each  constituent  entering  unit  volume  per  unit  time, 
either  due  to  convection  or  diffusion,  equals  the  mass  quantity  of  the  constituent  disappearing 
as  a  consequence  of  the  chemical  reaction. 

These  equations  are  as  follows 

ot>  vK,  1  V  (pK,e, 1 5] 

where  i\u  is  the  diffusion  velocity  of  species  i,  and  ir,  is  the  mass  production  rate  per  unit 
volume  of  species  i. 

We  consider  a  one  step  chemical  reaction  in  which  the  forward  reaction  is  dominant,  the' 
backward  reaction  being  negligible.  For  an  Arrhenius  type  reaction  with  second  order  che¬ 
mical  kinetics,  we  may  write 

Hyp  —  f;(/>/HT)  exp  |  K/I1T|-K,K.J  |(ia] 

where  H  is  the  activation  energy  of  the  reaction  and  h  is  the  frequency  factor.  Also  if  v  is  the 
stoichiometric  ratio  oxidizer-fuel 

ic,  »if,,  )/■,=  — (1  •  v)ny  [lib] 

We  shall  urn  relation  |(>a]  through  most  of  this  study.  The  extension  to  more  genera! 
reaction  rates  of  the  form 

ayo  -  — //(T,  />)expl —  K/HT)  KV  Ki  |7] 

is  easily  made. 

The  diffusion  velocities  depend  on  pressure,  temperature  and  species  concentration 
gradients.  1’sually  the  pressure  gradient  effect  on  diffusion  velocities  is  small  compared  to 
those  due  to  mass  fraction  gradients.  This  is  specially  true  when  mixing  takes  place-  in  thin 
mixing  regions  and  boundary  layers.  Thermal  diffusion  will  be  neglected  because,  diffusion 
velocities  due  to  gradients  of  temperature  are  generally  a  small  fraction  of  the  velocities  due 
to  concentration  gradients. 

If  t'  •  molecular  musses  of  die  species  are  approximately  equal  we  may  use  Kick’s  law 
for  the  determination  of  r,h. 

K,r„  DTK,  [«] 

where  D  is  an  average  diffusion  coefficient. 

If  the  concentration  of  one  of  the  species  is  small.  Kick's  law  is  valid  for  the  other  two 
species.  This  always  happens  in  diffusion  flames  where  oxidizer  concentration,  for  example, 
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is  very  small  in  the  reaction  zone,  or  in  the  fluid  side  of  the  flame.  Then  we  may  use  Kick’s 
law  for  fuel  and  oxidizer  with  the  diffusion  coefficients  determined  by  the  binary  mixtures: 
fuel-products  and  oxidizer-products  respectively.  In  this  study  we  will  use  a  single  average 
diffusion  coefficient  1). 

Inserting  [X]  into  [5]  we  obtain 

r  ■  V  K,  =  ( 1  'o)  V  •  (<(  I )  V  K,)  l-ir./p.  [!)] 

11.4.  Momentum  iqiuitiim. 

c  Vc=  (1/p)  V p  j-  (l/[<)  V  3.  [10] 

Where  a  is  the  stress  tensor 


We  neglect  the  diffusion  stress  tensor.  Gravity  forces  will  also  be  neglected  for  simplicity, 
although  they  can  only  lie  neglected  for  large  Fronde  numbers  and  this  is  not  always  the 
ease  in  diffusion  flames. 


:t  !l  dx. 


l  4-  f'’ 

\dxj  dx, 


:;)■ 


ll.f).  Kneniji  Ktfuuium. 


If  the  specific  heals  c,,,  of  the  species  are  assumed  to  be  ciptal  and  constant  the  energy 
eipiation  may  be  u  rillen 


VT 


'  « •  V  ft 

pc(, 


»'  1  </ 
P  C„ 


|H] 


where  t  :  Vr  is  the  Ihtyleigh  dissipation  function  and  r/  is  the  chemical  energy  that  a  combus¬ 
tible  mixture  containing  a  unit  mass  of  fuel  and  v  units  of  oxidizer  has  available  for  conversion 
into  thermal  energy  if—  It"  j- v/i!J  (1  f-v)/*JJ.  IV is  the  I’randtl  number  which  will  be  assumed 
constant.  Thermal  radiation  is  not  taken  into  account. 

Kipialions  [2|,  [4 1,  [10],  [11]  and  [!l[  (for  i— 1,2),  together  with  relations  [(>]  and  the 
functional  relations  between  the  transport  parameters  and  ft,  T,  and  K,,  constitute  the 
system  of  differential  equations  governing  the  structure  of  diffusion  flames. 

In  addition  we  must  include  the  appropriate  boundary  conditions. 

Without  losing  much  generality  we  can  state  as  boundary  conditions  for  the  tempeiature 
and  mas  fractions  that  they  be  constants  at  some  surfaces  or  zones  of  the  flow  field. 

For  example,  on  some  surface  or  region  at  the  fuel  side  of  the  flame  —  the  fuel  exit  — 


and 


K/~  K„  „  T  — Tr  j 
K,=  K .  T=T«  \ 


1 1 2] 


on  some  surface  or  region  lac  from  the  Maine  on  the  oxidizer  side  of  the  flame. 
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II. 6.  Dimensionless  Form  of  the  Equations. 

Let  us  introduce  the  following  non-dimensional  variables 
x  =  Lx *  „  i!=Un*  „ 


T=T„  +  '  T* 


!‘o|i 


Re== 


pop 

p0UL 

!Lo 


/)  =  p„l  I1;)* 

Se  =  11 
pl) 


[13] 


He  is  the  non-dimensional  Reynolds  number,  and  fie  is  the  fielunidt  number  that  will  be 
assumed  to  be  constant  and  equal  to  the  l’randtl  number. 

Subscripts  „  and  will  indicate  boundary  conditions  far  from  the  flame,  on  the  oxidizer 
and  fuel  side  of  the  flame  respectively. 

The  characteristic  length  Land  velocity  lT  are  some  overall  characteristic  magnitudes. 

In  terms  of  these  non-dimensional  variables  the  governing  equations  take  the  form. 


/»« 

P.  tT! 


o* = o*  / 1  |-|  ''  T* 

'  I  ('l>  ‘  t» 

V*.  ([/?*)=«) 


1 


1 

He 


V*-  o* 


r*  •  V'T* 


He  1 


(>*•  V*c*=  —  V*p* 

f>r 

r*K‘“|  Jiv  ] ^ v‘  1  |  oj(”o)’ 

I  ,!•  '•  rT->  - 1 I  ( 'T f  1 7 1  ■ r"' 1 1  </  it*)  I  c r 


>|14J 


Here 


where 


[15] 


'17  will  be  the  adiabatic  temperature  of  the  flame 


q  K,,.(T*  -I-  K,,  ) 
C|,  IvJn  i  vl\j,. 


as  we  shall  see  later.  Then  tc  is  a  characteristic  chemical  time,  such  that  (m, •/(>)*  will  ho  of 
order  unity  if  the  mass  fractious  are  not  small  and  the  temperature  is  close  to  the  adiabatic 
flame  temperature. 


[20] 
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Solution  [20]  is  independent  of  the  chemical  kinetics.  There  are  cases,  however,  in  which 
the  boundary  conditions  as  given  in  [12]  are  not  known  m  priori*  because  they  depend  on 
the  chemical  kinetics.  For  example,  in  the  ease  of  a  fuel  droplet  burning  in  an  oxidizing 
medium,  the  oxidizer  mass  fraction  at  the  droplet  surface  (which  is  zero  for  infinite  or  very 
largo  reaction  rates)  may  build  up  to  some  unknown  value  when  the  reaction  rate  be¬ 
comes  low. 


1 1 1.2.  Iliirke-Scltununni  Solution. 

For  the  study  of  diffusion  flames,  Burke  and  Schumann  introduced  the  assumption  'hat 
the  region  where  ivjr,  is  different  from  zero  is  infinitely  thin,  and  K,=0  on  one  side  of  the 
flame,  and  Ka  =  0  on  the  other. 

This  should  he  true  when  l./F  lr  is  very  large.  If  both  K,  and  Iv..  were  different  from 
zero  in  a  region  where  the  temperature  is  not  low  compared  with  T f.  then  (»•,/<>)*  would  be  of 
order  unity  and  the  term  [l./T  tr]  (wy/o)*  would  be  very  large  compared  with 

and  |  u-11,1.]^*.(p*V*Ki), 


that  are  of  order  unity. 

Also  if  in  system  |!4|  we  take  the  limit  l./l'/r->ao,  then  we  obtain  the  result  ( ny';, )*  =  ()• 
So  either  !<.,==<>  or  K._.  0,  and  system  1 14]  takes  the  following  form,  where  f=l,:i,  j='2 

on  the  fuel  side  of  the  flame  and  i= 2,3,  ./=-!  on  the  oxidizer  side  of  the  flame. 


::'r  k  ‘■•("IJdi 

V*.  (;,*»•*)  =  0 


r'.f'r' 


1  I  1 


He  ^ 


»•*•  V*T 


He  IV 


r.(n*H”i 


K, 


1  l'1 

He  q 
=0 


:  V*  r' 


[21] 


We  ran  use  with  system  [21]  the  same  boundary  conditions  [12]  (*)  of  system  [14],  if  we 
allow  for  discontinuities  in  the  mass  and  temperature  distributions  at  the  zero  thickness  flame. 

In  the  flame  the  captations  of  conservation  of  mass  and  energy  indicate  that.  1)  fuel  and 
oxidizer  diffuse  towards  the  flame  in  stoichiometric  proportions;  2)  that  the  heat  leaving  the 


(*l  K10  ;ui«l  K2»  must  In*  7.<»ro  if  I.  I*/,  ►  r. . 
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flame  duo  to  enmluction  equals  the  heat  released  by  the  reacting  species  when  reaching  the 
flame.  That  is 


<5k,  an, 

V  I  T~~  — - ■— 

du,  du., 

g  <?k,  ar  ^  a-r 

<■,,  dii,  dn,  dn.. 


[22] 


where  d  dn,  and  d/dn,  indicate  derivatives  normal  to  the  flame  surface  toward  the  fuel  and 
oxidizer  sides  respectively.  They  must  he  evaluated  at  the  flame. 

The  temperature,  mass  fractions,  and  therefore  the  density  and  velocity,  are  continuous 
functions  at  the  flame.  For  this  reason,  mass  and  heat  transport  by  convection  is  not  taken 
into  account  when  writing  the  conservation  equations  through  the  flame. 

lly  solving  the  system  of  equations  [21],  with  boundary  conditions  [12]  and  [22],  we  obtain 
the  Ihirke-Sehumann  velocity,  mass  fractions,  and  temperature  distributions. 

In  particular,  let 


K,  ■ —  Ks/v=/‘|(oi)  I 
K,  :  T*  —  /i,  (.<•)  ( 


[23] 


be  the  solutions  of  equations  [lfi]  and  [18]  (valid  only  for  V-/q^  1  and  (1/He)  F2/<7'4 !)• 
The  equation  of  the  flame  surface  is  obtained  by  writing 

K,  K.,  (I  /',  (r)  '■—(). 


Also  according  to  |20|  at  the  flame 


ami 


(K„  •  -T*)K.„, 
>K„.  I  K„ 

(K,,.  i-  t;»  k,„ 

•'Ku  i  K >„ 


[24] 


lly  writing  K,  —  0  on  the  fuel  side  of  the  flame  surface  and  K,=0  on  the  oxidizer  side, 
we  obtain  the  temperature  and  mass  fraction  distributions. 

K,  —f,  (•<•)  /  K,  =t)  j 

K,-=«  .for/',  n  K,  -  —  v  f,  (,r)  '•  for  /',  -  '  0.  [25] 

I  *  ~f<M)  —  f,  (•<•)  i  T*  f,  (x)  ) 


It  is  interesting  to  point  out  that  the  Ilurke-Kchumann  solution  satisfies  the  complete 
system  of  equations  [14]  and  also  its  boundary  conditions.  This  solution  is  not  the  correct 
one.  only  because  the  first  derivatives  of  the  temperature  and  mass  fraction  distributions 
have  discontinuous  first  derivatives  within  the  flow  field. 

Solutions  1 23]  of  equations  [Hi]  and  [IK]  are  modified  when  finite  values  of  L/Ff,  are  con- 
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siderod.  Tho  reason  for  these  modifications  is  that,  although  reaction  rates  do  not  appear 
explicitly  in  equations  [1(>]  and  [18],  the  variables  r  and  |i  that  appear  in  those  equations 
will  depend  on  the  reactions  rates. 

However,  wo  may  expect  that,  for  sufficiently  large  values  of  L/lI/„  tho  reaction  zone  (or 
region  where  u\  4“ 0)  will  ho  very  thin.  Hence  tho  Hurke-Schumann  solution  [25],  for  which 
the  reaction  zone  has  zero  thickness,  will  ho  a  very  good  approximation  in  tho  ease  of  large 
but  finite  I./'U lr.  This  will  ho  especially  true  outside  of  the  reaction  zone. 

Equations  [1(1]  and  [18],  in  particular,  should  remain  practically  unchanged.  This  is 
exactly  right  when  mixing  and  reaction  takes  place  in  constant  pressure  regions  and  boundary 
layers  if  o re  is  assumed  to  he  constant.  In  such  eases  equations  [1(5]  and  |18]  as  well  as  the 
boundary  conditions  (for  large  I, /l’  will  he  independent  of  the  reaction  rates.  The  same 
will  happen  then  to  their  solution. 

Summing  up:  If  the  reaction  rate  is  sufficiently  largo  the  reaction  zone  will  he  very  thin 
compared  with  any  other  important  length  (as  for  example,  the  width  of  the  mixing  region). 
Then,  in  order  to  obtain  the  velocity,  temperature  and  mass  fraction  distributions  outside  of 
the  reaction  zone,  i.  e.,  for  the  study  of  the  external  structure  of  the  diffusion  flame,  we  may 
use  the  assumption  of  infinite  reaction  rates. 


IV.  STitrrrt  itK  ok  tiik  khaction  zonk 
I V .  1 .  Till*  >  ('In' mini  I  llountlorii  l,tn/ir  . 

The  fact  that  in  the  limiting  ease  of  infinite  reaction  rales  the  thickness  of  the  reaction 
zone  is  zero,  and  that  (lie  first  derivatives  of  Kf  and  T  normal  to  the  flame  tire  discontinuous 
there,  suggests  that  for  large,  although  finite.  L/l’f,: 

a)  The  thickness  of  the  reaction  zone  will  he  small, 
h)  The  diffusion  terms 


1  d 

!'  ^  h,  \ 

and 

LA, 

1  |i  <*T\ 

;<  dn 

IT  dn  I 

o  dii 

i  IT  dn 1 

will  balance  the  chemical  production  terms  irp'o,  these  terms  being  very  large  compared 
with  till  the  other  terms  of  the  equations.  (Here  d  dn  indicates  differentiation  normal  to  tho 
flame). 

In  other  words,  for  large  values  of  the  chemical  reaction  rates  the  reaction  zone  will  he 
a  very  thin  region  or  chemical  boundary  layer  .  There,  due  to  the  rapidly  varying  gradients 
of  temperature  and  mass  fractions  normal  to  the  flame,  mass  diffusion  and  heat  conduction 
normal  to  the  flame  constitute  the  only  transport  mechanism  required  to  balance  the  chemical 
production  terms.  Transport  by  diffusion  or  conduction  in  other  directions  or  convection 
may  be  neglected  within  the  reaction  zone. 

In  order  to  show  this,  let  us  assume  that  we  know  the  Itiirke-Sehumann  solution  [25], 
lienee,  we  know  the  flame  surface  location  for  infinite  I.  l'/r.  and  therefore  the  approximate 
location  of  the  flame  region  for  large  I,  f/r. 

For  simplicity  we  w  ill  limit  ourselves  to  the  two-dimensional  ease.  Tho  results,  however, 
are  completely  general.  We  shall  write  the  equations  of  motion,  and  mass  and  energy 
conservation  equations  in  a  curvilinear  system  of  coordinates.  In  this  system,  see  Fig.  1, 
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x  will  lie  the  distance  measured  along  the  flaino  surface,  as  determined  by  the  Hurke-Schu- 
mnnn  solution.  The  distance  normal  to  this  surface  will  be  indicated  by  g;  u  and  r  will  be 
the  velocity  components  in  the  x  and  g  directions;  1/K  is  the  radius  of  curvature  of  the  flame 
at  point  x.  Limiting  ourselves  to  a  region  where  K //  is  small  compared  with  1  the  line 
element  has  components  (1  4-  Kg)(/.r  and  rfg.  and  the  eipiations  are  as  follows: 

('iiiiliuiiil/l  for  the  mixture 
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♦r=M2(*J,  4,d  4,  2'Hu^  i- 

and 

I  d  ii  K  r 
I  K  //  d.r  t  !  Ivy/ 

d  r 

**  d.n 

1  dr  dn  K  n 
I  K  //  d.r  d/i  I  :  K  I/ 

Now  lot  1,  ami  r  lio  overall  characteristic  leiifith  and  velocity.  I.ot  ami  6m  bo  (lie 
thicknesses  of  tlio  roaotion  /otto  ;m< I  of  (lie  mixing  region  respectively.  Wo  could  show  tluit. 
for  low  lleynolds  innnhor.  I.,  while  iM  |  l)„l.  lT  for  larfie  lioynolds  mtnihors. 

The  mass  fractions  just  at  the  outer  odjio  of  the  reaction  zone  will  ho  of  order  The 

satin'  order  of  niaonitudc  w  ill  ho  valid  inside  of  the  reaction  zone. 

hot  us  now  introduce  the  nnii-dinicnsionnl  variables  a'/l.,  ///or,  (trap)  ,)s,  by 

slrelehino  the  coordinates,  mass  fractions  and  reaction  rates,  so  as  to  make  (he  non-dimension¬ 
al  factors  and  terms  of  order  unity  within  the  reaction  zone.  For  the  remaining  variables 
we  may  use  the  same  non-dimensional  variables  used  in  section  II. 

We  will  write  the  eoverniii"'  equations  in  terms  of  these  non  dimensional  variables. 
Now.  if  the  terms  account  ini:  lor  the  conduction  and  diffusion  normal  to  the  flame  are  o'oinje 
to  he  of  the  order  of  the  chemical  production  term,  o;' 

If  we  now  lake,  in  tin-  ouerjiy  and  diffusion  equations,  the  limit  <<, ,  Sl(,  •(),  most  of  (he  terms 
in  these  equations  drop  out.  We  are  left  with  the  following  differential  equations: 
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:•  d  K,  \ 

;<  d  /i  \ 
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:■  <5T  t 
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dn  l 

I’r  d/i  I 

that  we  have  written  in  dimensional  form. 

l-Yoin  the  inomenlinu  equation  we  deduce  that  the  variations  of  pressure  across  the 
reaction  zone  are  of  order  •),  l„ 

lienee  we  may  assume,  when  writing;  the  equation  of  state,  that  the  pressure  is  constant 
across  the  reaction  region  and  equal  to  the  value  obtained  at  the  flame  with  the  llurke  Schu¬ 
mann  assumption. 

Therefore  hi  addition  fo  the  above  equations  for  K , .  K...  and  T  we  have  the  equation 

/»(./■)  ;,TI!M 

where  />|.rt  is  a  known  function  of 
Also 

ir..  ir.  Ii  ■  |  /<(.r|  UT|  exp  . —  KI(T:K,K,  |'2!t] 

or  we  may  use  tin-  ueneral  expression  |7|  for  u,  as  is  done  in  appendix  A. 
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N'o  derivatives  with  respeet  to  ,c  appear  in  system  [2<i]  to  [2!)].  Therefore  these  equations 
may  he  solved  as  ordinary  differential  equations  in  whieh  a-  stands  as  a  parameter. 

As  boundary  renditions  we  will  write  that  when  //-»  +  oc  the  temperature  and  mass 
fraelion  distributions  eoineide  with  those  obtained  by  assuming  the  reaetion  rate  to  he 
infinite. 

IV.2.  The  Solution  of  llic  Chemical  Hoinnlon/  Lui/cr  h'linalionx. 

r 

l!y  iatrodueint;  the  new  variable  //,=  /  (--/(>„)<///,  if  we  assume  that  p ji  =  o„ p0,  equa- 

•-  o 

tions  [27 )  and  [28]  may  he  written 


<1-  K;  1  II'. 

•ini  1  ;> 

[30] 

d-T  l  ir, 

<lifi  l',,  [• 

l»l] 

From  equations  |:t()|  if  we  take  into  account  that  ir.,—  vw,  we  e(q 

K,  K;  »~A|//|  :  It,. 

l»2]  1 

Similarly,  from  [30|  and  |:tl|  we  obtain 
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‘  J 

K,  T*  A.,//,  ■  It.. 

j 

K,v  T*  A ■  1 

l»4]  ! 

delations  [:I2|  to  [3-1 1  are  independent  <d'  the  ehemieal  reaction  rates.  However,  they  are 
only  valid  within  the  reaetion  /.one.  The  constants  A,.  I!,  must  hi*  chosen  so  that  these 
relations  eoineide  with  the  similar  relations  obtained  from  the  Iturke-Sehunianu  solution,  at 
least  for  low  values  of  //,.  Then,  relations  [33]  and  |:M]  may  lie  written. 
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Now,  by  using  expression  [fin]  for  the  fuel  production  rate,  equation  [31]  takes  the 
following  form: 

rf-'fi  --'I y  P  i  \  o„  *  o„  i  — e  (|  v  //,  1 1  i  i  v  :  !h  i 

—  . rr  •  ■  •  i  1 1  “ 1  +  k„  d  1 1  hr  ~  id  K  I  W 


'hi  i 


1U,  p„  0„  !-0  /  1  I  ff,  fl„  jfi 


where  Dn~c,,T„lq  and  ()ll  =  Kc,,:Hq  is  the  non-dimensional  activation  energy  of  the  reaction. 

Ity  solving  equation  |3(>]  with  the  boundary  conditions  K, =0  for  y/,  — » oo  and  K,=0  for 
i/,  — ►  — oc,  we  obtain  the  temperature  distribution  within  the  reaction  zone. 

If  the  reaction  rate  is  sufficiently  large,  the  temperature  will  not  deviate  appreciably, 
within  the  reaction  zone,  from  its  limiting  value  (when  le— *())#==  1  at  //,  0. 

Then  a  good  approximate  solution  of  [3(i]  may  he  obtained  by  substituting  the  factor 

I  .  n  l  — «  | 

exp  ■  t 

«!„  !  1)  1  )  14-  U„  «„  4-  0  i 

by  its  value  at  //,=<>.  Let  « 0,. (.»*)  for  //=().  Then  we  approximate  equation  [:i(i|  by 


Kqualion  |:i!l|  was  solved  numerically  (15).  Its  solution 


|38a] 


is  represented  in  Fig.  with  a  solid  line.  In  particular  ,1(o)  O.KIili.  An  approximate 
solulion  is  presented  in  Appeiutix  A. 
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This  we 

may  write  in  the  form 
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where 


A. = 4  v  T,«  *l  ''L  - = 4  v  T/W  —  ‘ )" 

'  />„  ! >»  tr  )>(,  w-tr 


A  second  approximation  for  valid  for  values  of  Or>  O.S  is 
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IV.:l.  Disnts-siuii  of  the  solution. 

liclalion  |4f>]  shows  that  deviations  of  (hr  temporal uiv  from  its  asymptotic  llurke-Srhu- 
inann  value  0,-—l  depend  on  the  parameter  A„.  This  parameter  incorporates  both  (ho 
chcmieal  kinetic  parameters  (through  the  value  of  the  elmraeleristie  chemical  time  I.)  and  the 
fluid  mechanical  parameters  (through  the  mixing  lime  ');■■,  (.<•)  1  )„).  A„~ 

The  obvious  result  is  that  the  deviations  of  the  temperature  and  mass  fractions  from 
their  limiting  asymptotic  values  increase  w'th  decreasing  values  of  the  ratio  tjic.  Now  tm  is 
inversely  proportional  to  the  fuel  rate  of  supply  to  tin*  flame  in.  lienee  we  deduce  that  by 
increasing  the  fuel  rate  of  supply,  the  flame  temperature  will  decrease. 

The  initial  fuel  and  oxidizer  mass  fractions  influence  the  results  through  the  factor  T/*!’ 
that  appears  in  |  l  !|.  !!y  diluting  the  fuel  or  the  oxidizer  we  gel  larger  deviations  from  the 

liitrkc  Schumann  result. 

Itclnlion  |  Hi]  indicates  that  the  thickness  of  the  reaction  region  decreases  with  increasing 
values  of  the  mass  rate  of  fuel  supply. 

Kor  this  chemical  boundary  layer  scheme  to  he  valid  i,.  must  lie  small  compared  to 
unity,  hence  the  same  criterion  may  he  used  for  the  validity  of  the  Ilnrke-Schumaim  assump¬ 
tion  and  solution,  and  for  tin-  validity  of  the  chemical  boundary  layer  solution. 

Obviously,  the  most  important  characteristics  of  the  chemical  boundary  layer  solution 
and  those  of  the  liurke-Scluununn  solution  coincide.  For  example,  flame  location  and  the 
■  pianlity  of  fuel  burning  per  unit  flame  area  are  the  same  in  both  solutions.  The  chemical 
boundary  layer  solution,  however,  gives  a  finite  thickness  for  the  reaction  zone,  and  a  small 
correction  to  the  temperature  and  mass  fraction  distributions  This  correction  can  he 
evaluated  very  easily  and  accurately  in  terms  of  the  chemical  kinetic  parameters  of  the 
reaction.  So.  the  chemical  boundary  layer  solution  can  be  of  help  for  the  study  of  chemical 
kinetics  by  means  of  experimentation  in  diffusion  flames  (I2|.  (28). 

The  parameter  A...  that  measures  the  deviations  from  the  infinite  reaction  rate  solution, 
max  be  used  for  the  determination  of  an  extinction  criterion.  This  is  supported  by  the 
follow  ing  reasons: 

a)  I  tile  to  the  high  values  of  the  activation  energy  of  many  of  the  chemical  reactions, 
the  chemical  production  term  is  very  sensitive  to  temperature  variations.  This  accounts  for 
the  hurt  that  llatne  extinction  occurs  in  a  rather  sharply  defined  way.  This  may  also  he  due 
to  tin-  existence  of  some  ignition  temperature  for  the  reaction. 
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l>)  The  concept  of  an  overall  reaction  rate  is  not  in  general  valid  through  a  large  tem¬ 
perature  range.  Then,  the  idea  of  solving  the  complete  exact,  equations,  for  obtaining  an 
« exact-  extinction  criterion,  loses  part  of  its  interest  if  use  lias  to  be  made  of  some  assumed 
overall  reaction  rate  expression  throughout  the  whole  temperature  range. 

c)  The  chemical  boundary  layer  solution,  although  it  cannot  explain  extinction  except 
in  a  qualitative  way.  can  provide  a  criterion  for  extinction  not  to  occur  if  the  overall  reaction 
rate  is  known  to  he  valid  in  a  given  temperature  range. 

Whenever  the  parameter  A„  is  sufficiently  low  as  to  make  T,..._.  O.NT,  (and  this  occurs  for, 
roughly,  A„  tit)),  the  thickness  of  the  reaction  region  begins  to  be  comparable  with  the 
thickness  of  the  mixing  region.  Then  the  rate  of  fuel  consumption  begins  to  diminish  and 
hence  T,  will  begin  to  decrease  even  faster  with  decreasing  A„. 

Therefore,  A„  SO  may  be  used  as  an  approximate  extinction  criterion  as  well  as  a 
criterion  for  the  validity  of  the  llurke-Schumann  solution. 

We  have  set'll  that  A„  -  pj  I )„/«»-’  tr.  Now,  this  same  parameter  appears  in  the  theory  of 
promixed  laminar  flames.  There  it  takes  a  value  of  the  order  of  KM),  that  depends  on  the 
initial  mass  fractions  and  energy  of  activation  of  the  reaction  (24),  when  ii«  is  substituted  by 
(In'  fuel  consumption  rate  (1  per  unit  area.  Therefore  an  ■approximate*  relation  may  be 
established  (!)),  (10)  between  the  value  of  m  at  extiuetiou  (maximum  flame  strength)  and  ii: 

^  (i. 


V.  Tun  hxtkknai.  sTiit’crruK  os  i.aminak  on  st  siox  si.amks 

V.l.  /.ioi/c  Hi'i/Holtls  Number  (Vise,  hiviticid  Kquiitioux. 

The  chemical  boundary  layer  equations,  governing  the  temperature  and  mass  fraction 
distributions  within  the  reaction  zone,  have  been  solved  in  till'  most  general  case.  However, 
for  the  detailed  evaluation  of  the  solution  we  must  know  some  parameters  appearing  there. 
They  include  the  mass  rate  iii  of  fuel  consumption  per  unit  flame  surface,  and  the  ideal  flame 
location. 

In  order  to  evaluate  these  parameters  as  well  as  the  mass  fraction  and  temperature 
distributions  outside  of  the  reaction  zone,  the  llurke-Schumann  solution  must  be  obtained 
first. 

As  mentioned  in  the  introduction  this  solution  has  been  obtained  in  some  particular  cast's. 
1'nfortnnately.  even  when  using  the  Iturke-Schumanu  assumption  of  infinite  reaction  rates, 
the  resulting  system  of  equations  [21]  anil  boundary  conditions  [12]  and  |22]  is  so  complicated 
that  only  a  few  approximate  solutions  exist. 

Marble  and  Adamson  (17)  have  pointed  out  that  a  number  of  important  combustion 
problems  may  lie  investigated  analytically  with  tin1  help  of  boundarv  layer  approximations. 
Most  of  the  solutions  so  far  obtained  make  use  of  these  approximations.  These  may  be  used 
whenever  the  Kovnolds  number,  based  on  some  overall  dimension  of  the  flow  field,  is 
sufficiently  large. 

\\  e  w  ill  show  here  that,  by  using  some  additional  assumptions,  a  fairly  simple  solution  of 
the  Ihirkc  Schumann  mixing  problem  is  obtained 


1 
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If  in  system  [21 J  we  take  (he  limit  He  *  *,  wo  obtain  the  following  system  of  differential 
equations,  that  \v<>  slmll  write  in  dimensional  form: 

j»=PTlf/.M 

t.(;,r)  =  0 

r  -  V  r  — - —  V  p 

a 

r  ■  r  K,==  <i 
•  t  > 

r  •  TT=  »'-r« 

?*V 

The  lioundary  conditions  [22]  cannot  he  reiained  because,  in  the  process  of  tailing  the 
limit  He  -~x,  we  dropped  the  higher  derivatives. in  the  equations.  j 

However  the  lioundary  conditions  [12]  emi  he  satisfied  if  we  allow  for  the  existence  ef 
discontinuities  in  temperature,  mass  fractions,  density,  and  velocity  at  some  stream  surface. 

The  position  of  this  surface  is  determined  l>y  requiring  that  all  the  boundary  conditions  |!2j 
lie  satisfied. 

if  we  consider  only  low  Much  number  flows,  then:  | 

I)  Tlit'  density,  temperature  ami  mass  fractions  will  be  constant,  although  possibly  with  ‘ 

different  values,  on  each  side  of  the  discontinuity  surface,  S <><«  pig,  t.  I 

21  Kqiiations  [17 J  reduce  to  the  system  i 

V  •»■-<>  j 

•  '  *  I 

eye  — VP  \ 

b  1 

and  tangential  discontinuities  of  r  are  allowed  for  at  sonic  surface,  so  ns  to  satisfy  (ho  j 

required  lioundary  conditions  on  r.  The  pressure  mast,  of  course,  lie  continuous  at  tlio 

surface.  ! 

As  an  example.  Hie  solution  of  this  problem  for  the  low  speed  source  flow  is  presented  in 
Appendix 

Y.2,  Mifithj  La/ fry  lu/Hatiwis. 

For  large  but  finite  He.  the  ideal  discontinuity  surface  is  substituted  by  a  thin  mixing 
layer  with  the  same  approximate  location.  Although  the  discontinuities  in  the  temperature, 
mass  fractions,  etc.,  no  longer  exist,  the  derivatives  of  these  variables  normal  to  the  mixing 
layer  will  lie  very  large  eompared  to  the  derivatives  in  tite  surface  direction: 

In  order  t*>  study  the  structure  of  the  mixing  region,  we  will  write  the  system  (21  j  in 
boundary  lay  er  coordinates.  Then  we  can  obtain  Hie  mixing  layer  equations  by  using  tho 
well  known  boundary  layer  limiting  process. 


i; 


j 

[48] 
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Limiting  ourselves  to  the  two-dimensional  or  axially-synnnetrie  low  Mach  number  flow 
cases,  we  will  write  these  equations  in  the  form  given  by  Lees  (25). 


v-  /  >>/?„)  <n- 

'=  /  W'dM-b/X 

J  n 


m 


where  A'=0  for  two-dimensional  flows  and  A-=l  for  axially-symmetric  flows;  u„ (,<■)  is  the 
velocity  at  the  oxidizer  side,  just  outside  of  the  mixing  layer. 

By  introducing  the  stream  function  'I-  such  that 

d'l-  dT 

<•  a  i*  ~=  ,  Y  !*  i'  —  •  [50] 

The  continuity  e<|uation  is  automatically  satisfied.  Let 

M'(S,rl)  =  l  2S/,(r„5).  Then  «/«„  fir,,*)  [51] 


whore  the  primes  denote  differentiation  with  respect  to  r,.  We  will  assume  fn  =  p(1|i0.  Then 
the  mixing  layer  equations  take  the  form 


r 


i ITT,,  ,™=u„.5  j 

K”  i  I'r/'Kj  -  =  t.  i.d.;  „  Ky  0  I 

T"  •  l’r/T^t.i.d.S  ' 


|52] 


Where  j  -.=  2  on  the  fuel  side  of  the  flame,  and  i=2, 0;  j=  1  on  the  other  side. 

Tin-  l.i.d.;'s  in  the  right  hand  side  of  the  equations  indicate  terms  involving  derivatives  with 
respect  to  s. 

As  boundary  conditions  we  may  write 


Where  r(/(-l  gives  the  ideal  flame  position. 
In  addition 

I'  1  ..  K..  =  K.,„ 

f--  ",  “=  K„ 


T  —  T„  for  r,  .  oc  / 
T  =  T,  for  r,  -* — oc  \ 


[58b] 


l-'or  the  solution  of  the  above  mixing  problem  a  third  boundary  condition  for  f  is 
required.  This  should  be  derived  from  the  compatibility  condition  of  the  higher  order 
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approximation  (2G).  However,  for  our  purposes  we  may  write  as  third  boundary  condition 
/■(0,  !)=<),  because  the  only  effect  of  changing  the  value  of  f( 0,9  is  a  displacement  of  the 
mixing  layer  in  the  Y  direction. 


V.!i.  Local  Similarity  Approximation. 

As  the  boundary  conditions  are  independent  of  the  functions  /',  K,  and  T  would  ho 
functions  only  of  r,,  and  similarity  would  exist,  if  the  presure  gradient  parameter  (2  £/«,>)  d  «„/</  s 
were  constant. 

This  occurs  at  the  stagnation  point  where  the  parameter  has  the  values  t/2,  for  the 
axially-symmetrie  ease,  and  t  for  the  two-dimensional  one. 

Similarity  also  exists  in  the  constant  pressure  case  corresponding  to  the  mixing  of  two 
parallel  streams.  In  this  ease  the  pressure  gradient  term  is  obviously  zero.  The  local 
similarity  approximation  limy  he  used  when  (2 ••/'<»)  </«„;</;  is  a  slowly  varying  function  of  5. 
Then  wo  may  neglect  the  t.i.d.-'s  in  the  right  hand  side  of  equation  |r>2|.  The  resulting 
system  of  differential  equations  may  then  he  integrated  as  a  system  of  ordinary  differential 
equations. 

It  is  interesting  to  point  out  that  the  factor  |T/T„  - (/")'-']  in  the  pressure  gradient  term 
approaches  zero  it t  both  edges  of  the  mixing  layer.  The  neglect  of  the  pressure  gradient 
term  is  similar  to  the  neglect  of  free  convection  in  diffusion  flames. 

The  pressure  gradient  term  is  neglected,  without  much  justification,  by  Spalding  when 
studying  the  opposed  jet  diffusion  flame  (12). 

If  our  main  purpose  is  assessing  the  effects  of  chemical  kinetics  in  diffusion  f Iannis,  only 
an  approximate  analytical  solution  of  the  equations  is  required. 

This  may  he  obtained  easily  if  we  neglect  the  pressure  gradient  term  in  the  momentum 
equation  [52].  Then  it  reduces  to  the  Hlassius  equation 

r  ■  rr~n  nj 

with  the  boundary  conditions 

/'(«)  0.  /''(*)=!.  /'(  *)-•).-«,  it,,. 


An  approximate  analytical  solution  of  this  equation  is  gi\cu  in  Appendix  it. 
approximation  is 


/  1/  it,, 


The  first 


lr,r>I 


From  equations  [52j  we  deduce  the  following  system  of  equations 
(K,  K.vf  I’r/'tK,  -  K,  v|"  =  0 

(K,  T*f  IT/tK,  T*l  11 


m 

[57] 
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Taking  into  account  the  lioundary  conditions  [53  h]  we  obtain,  if  I’r  1, 


( K,  \  —  K , )  —  ( K  ,„  v  —  K, 


lv  |  ,.  K.ii.  ■ 

T*  i  K,  1 
Tf  ;  K„ 


1  —  erf 


.(MU 

I",1-) 


m 


And  the  boundary  conditions  at  r„  arc  satisfied  identically. 
The  flame  surface  is  located  at  the  point  r,=rlf  where 


,..f  /  1 1  ■  |  _  I'l  v  Kjo/v 

\  2  ~  Ktr  i  K.,„  v 


and 


in  —  I  o  I ) 


d  K, 


di/  | 


(7,'L-T.r  "it;  " 


p.-;  i  i$K<\  =  psi>» 


i  i  * 


[5-1] 


m 


So 


1  i*  ..  ,  1 1  i  >•  i  1  I  >■ 

i„, ,  K ,  |25  (K"  K"'')  2|=  ,>X1,i  4  V 


[til] 


The  parameter  A„.  measuring  (lie  deviations  from  the  lhirke-Srhimiann  solution,  can  now 
lie  evaluated.  Also  the  temperature  and  mass  fraction  distributions  outside  of  the  reaction 
•/.one  may  he  obtained  from  |57]  and  [58]  by  putting  K,= 0  for  r, rtl  and  Kj  =  ()  for  r,-  rlf. 

In  the  particular  ease  ).=  1,  we  get  /'=!,  and  tin*  equations  may  be  solved  even  for 
I’r  i  1.  We  obtain 


K,  K .  v  =  K , ,  <  K , ,  -  K,,„  v) 


1  !  erf 


I’r 


T  •  K,--:|T?  •  K , , ) 


1  erf 


i’r 
*>  r‘ 


iii 


S I  >„  a,.  il  .. 

I  I'm 


I'  |  |  I'  MI  ' 

Ktl  K,„* 
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Temperature  at  the  Burke-Schumann  flame  surface  as  a  function  of  the  distance  to  the  initial  mixing-  point. 


V.4. 


I’tirlicular  (''ixrx. 


In  particular  lot  us  consider: 
a)  Axially-symmotric  stagnation  point 


"i  ___  I  I',  1>',. 

I  I',. 

And  accordin';  to  1 44 1 


1>)  Mixing  ol  two  parallel  streams 

l>„ 


m 


It  „.r 


•\„  tvT.*->  ; 


K; 


Tlioroforo 


1  K._,„  v)  j  ( '4  _!  l)  exp  j 

1  i  r.  ( . 

4  T|'  ( 

[63J 

An  U„l 

(64] 

..  .  .  1  1  i  \ 

i  l\..„;v)  exp  - 

2 1  2  -  1 

1  1  >•  ,1 

4  Vj 

[65] 

(K,,  •  K.„,  •/)■-' exp! 

!  >•  / 

4  V  j  ■ 

[66] 

A.  rV„lr. 
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In  reference  ( 15)  a  numerical  integration  of  the  equations  for  the  mixing  ami  simultaneous 
chemical  reaction  of  two  parallel  streams  of  fuel  and  oxidizer,  moving  with  the  same  velocity, 
was  carried  out  in  order  to  compare  with  the  chemical  boundary  layer  solution.  The  results 
are  shown  in  Kies.  4.  5  and  li. 


VI.  ItKSf.MK 

Kxpcrimenls  (2)  and  the  success  id'  the  existing  theories  on  laminar  diffusion  f  lnmos  have 
clearly  shown  that,  in  those  cases  where  a  stable  laminar  diffusion  flame  has  been  obtained, 
the  Ifurke  Sehumann  assumption  (infinitely  fast  reaction  rate)  applies.  However,  the  Bur¬ 
ke  Schumann  solution  is  independent  of  chemical  kinetics,  and  does  not  give  any  criterion 
either  for  flame  extinction  or  for  the  validity  of  the  solution. 

The  fact  that  in  this  solution  the  flame  thickness  is  zero  suggests  that  in  practical  eases 
the  reaction  zone  must  lie  of  negligible  thickness,  making  it  possible  to  obtain  a  solution  of 
tile  boundary  layer  type  including  the  effects  of  chemical  kinetics. 

At  each  side  of  the  extremely  thin  reaction  zone,  chemical  reaction  effects  are  neglected 
as  compared  w  ith  convection,  conduction,  and  diffusion  effects.  The  reaction  zone  reduces 
to  a  flame  front  of  negligible  thickness,  which  acts  as  a  sink  for  the  reactants  and  as  a  source 
for  the  heat  and  products  evolved  in  the  chemical  reaction.  The  location  of  the  flame  front, 
rate  of  burning,  and  temperature  and  concentration  distributions  in  the  exterior  of  the 
reaction  zone  are  determined  by  using  the  Burke  Schumann  assumption. 

In  order  to  anal/.ye  the  structure  of  the  burning  zone  we  nitty  neglect  in  it  convection 
elici  ts  as  compared  with  conduction,  diffusion  and  chemical  reaction  effects.  The  equations 


governing  lliis  chemical  boundary  layer  are  ordinary  differential  equations  with  boundary 
conditions  determined  by  the  1  lurk  e-Selm  inann  solution.  Temperatures  there  are  elose  to 
the  adiabatic  flame  temperature,  and  then  an  overall  kinetic  scheme  applies. 

The  criteria  for  extinction  of  the  flame  and  for  the  validity  of  llurke-Sehumann  as¬ 
sumption  approximately  coincide,  and  may  be  obtained  by  solving,  once  the  llurke-Sehumann 
solution  is  known,  the  chemical  boundary  layer  equations.  This  solution  also  provides  the 
temperature  and  concentration  distributions  in  the  reaction  zone. 

If  our  main  purpose  is  the  evaluation  of  the  chemical  kinetic  effects  in  diffusion  flames 
an  approximate  solution  of  the  llurke-Sehumann  mixing  problem  will  be  sufficient.  This  we 
may  obtain  easily  for  large  Itoynolds  numbers  by  using  well  known  boundary  layer  ap¬ 
proximations. 
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APPENDIX  A 


In  this  appendix  an  integral  method  for  the  approximate  solution  of  equation  [81]  is 
presented.  Wo  will  use  the  same  approximations  as  those  used  in  obtaining  [37]  front  [36]. 
Wo  will  consider  the  more  general  mass  production  rate  expression  given  by  |7].  Then,  if 
relations  [35]  are  taken  into  account,  hv  taking  [7]  into  [31],  we  obtain 
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Here,  however,  we  will  not  choose  T,.(.r)  as  the  temperature  tit  //,=<),  hut.  as  the  tempera¬ 
ture  at  a  point  that  we  will  determine  later  on, 

Hv  introducing  the  variables 
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eiptation  ].\.i]  transforms  to 
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[A.5] 


That  must  be  solved  with  the  boundary  conditions 

1  o  tor  *  ‘  S 


In  order  to  solve  eipialion  jA.5j  in  an  approximate  way.  we  will  use  the  following 
integral  method. 

If-,  is  the  point  where  "T  reaches  its  maximum  value,  the  right  hand  side  of  equation  [A.5] 
may  be  written  in  the  form 

Where 

>.  tl-,i  ami  In:) 
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Neglecting  higher  order  terms  in  the  expansion  of  /*(')  in  powers  of  (:  —  r,).  equation  [A.B] 
may  lie  written  in  the  approximate  form 
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Tile  constants  ft,,  r,  and  k  will  he  determined  as  follows: 

a)  The  solution  jt(r)  of  [A.B]  should  satisfy  the  following  houndary  conditions  for  ft(«r); 
namely  ft'  —  *  t  1  lor  -  »  +  *. 

Then 
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and  integrating  [A.B]  between  c,  and  o c ,  we  get  fi',  =  0  (*). 
1  >)  If  by  means  of  the  relation 
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we  define  ft--fi(r).  then  we  will  choose  the  parameters  ft,,  z,  and  k  so  that  the  following  rela¬ 
tions  hold: 

!*'(*■■>  =  <>  anti  ft"(zt)  (ft,  —  *,)"  (Pi  +  «i)''-  [  A.9] 


h'rolli  relation  |A.8|  we  deduce 
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In  particular  for  r  — .r,.  we  obtain 
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(*)  H  in  addition  tiso  vtould  In*  made  of  tin*  boundary  condition  £  ♦  z  =  0  for 
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?  =  ?,  ■  --•'i-ft--  1  J  1  i'X|i< 

J  k  ~ 


w  hicb  is  also  a  fairly  t^tintl  approximate  -  dutimi  of  t»  [nation 
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From  (A. 10]  wo  deduce 
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Then  from  ]A.!t].  (A. 11]  and  (A. 12]  wo  get  the  relation 
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that  has  the  solution,  if  l/(«-|  b  1  (•)=,«, 
.  (i  j  ft  ]  2 


I  Je(| 


1  + 


*(«  I  b) 
2a  b 


■ir 


Now  c,  and  k  may  ho  doiormined  in  terms  of  3t  by  moans  of  the  relations 
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As  an  approximate  solution  of  (filiation  |A.1]  wo  may  use  either  the  function 

?(~>  ?,  ‘  <-  r,)orf/.  (c  r,)  ( |  ^  A)“ 1  ;  1  oxp(  k- („-  ^,1'-’)!  [A. 1(1] 

ohlainod  by  intopratinp  otpiation  |A.(i|.  or  the  more  approximate,  but  also  more  difficult  to 
evaluate,  function  3(")  solution  of  the  algebraic  oipiation  |A.K|.  In  both  oases  the  parumetors 
3i. and  /.■  are  as  fjivon  by  | A .  1 4 1  and  |A.15). 

in  order  to  compare  the  above  approximate  solutions  of  eipmtion  [A. 5]  with  the  exact 
numerical  solution,  we  consider  the  particular  ease  a  1.  It—  I.  Then 


3, --=0.87. 


=  I).  /.•  (U17 


and 


3U1---0.87  i’erf  0.(57  r  0.K4 1 1  exp  |II.(I7;)!] 
3(~'|  |r-  •  O.Ttioxp  (11.(17 


[A. 17] 

IA.I8I 
|  A.  1 51] 


The  approximate  solution  3(~~)  lias  been  plotted  with  a  dashed  line  in  Kip.  2  to  compare 
with  the  exact  solution  drawn  with  a  solid  lino. 

(loin;;  hack  to  the  penoral  ease  in  which  n  and  b  are  not  necessarily  one,  let  us  choose 
the  temperature  at  the  point  r,  as  the  temperature  used  in  approximating  equation  |H1). 

Then  from  relation  [A.:l|  by  makinp  r  r,.  we  pet 
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This  is  an  alpohraio  oipiation  for  b  . 
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Again  a  first  approximation  to  the  solution  is  obtained  by  making  %  1  when  evaluating 

y.(x,f>r).  Then  relation  (A. 'JO]  transforms  to 
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A  second  approximation,  similar  to  4?>,  may  he  obtained  for  f)r  without  major  difficulties. 
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In  order  to  approximately  solve  the  lllassius  equation 

r '  ir  » 

with  the  boundary  conditions 

/"<*>  i-  /■(--  x  )=->..  /■(<>)= o, 

following  Meksyts  (27)  let  us  write  equation  [H.l]  in  the  form 
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where  ni  =  f{o)  and  ii^  —  fio)  are  assumed  to  be  known. 
Then  by  integrating  [11.2]  twice  we  get 
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As  the  shear  stress  decreases  very  rapidly  with  increasing  values  of  [ v,  j,  we  may  expect 
that  most  or  the  contribution  to  the  integral  [H.:l]  comes  from  those  low  values  of  jr;|,  for 
which  f  can  be  well  approximated  using  only  the  first  terms  in  the  power  series  expan¬ 
sion  [H.2|. 

If  we  write 
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Then  the  integral  [B.3]  may  be  evaluated  and  we  obtain 
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If  wo  now  make  ►  +  x  ,  we  obtain 
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The  sories  in  [H.(>]  may  be  expected  to  converge  very  rapidly.  Keeping  only  the  first  two 
terms  of  the  series  expansions  appearing  within  brackets  in  [1US).  we  get. 


and 


that  may  be  approximated  with  less  than  4"  „  error  by 
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Now  by  means  of  relations  |li..X|  the  solution  |l!.'i|  may  be  written  tip  to  the  second  ap¬ 
proximation.  A  first  approximation  is 
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APPENDIX  (' 


Low  Si*kki)  Horitoii  l't.ow. 


Thi'  laminar  diffusion  flame  produced  wIhti  a  point  source  of  fuel  is  immersed  in  an 
oxidizer  stream  has  been  qualitatively  dismissed  by  Penner  (28). 

Let  the  density  and  temperature  of  the  fuel  leaving:  the  source  he  equal  to  the  density 
and  temperature  of  the  oxidizer  stream.  If.M  is  the  mass  rate  of  supply  of  fuel  due  to  the 
source,  the  velocity  distribution  is  given  by 


;  v(.v*v  4  *  ) 
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where  i ,  -  j  .r-  •  .cri  •  and  the  .r..  axis  is  in  the  free,  oxidizer,  stream  direction.  The  source 
is  iocated  at  the  origin  and  l’r  is  the  free  stream  velocity. 

The  stream  surface  that,  separates  fuel  and  oxidizer  is  easy  to  calculate,  and  that  has 
been  done  elsewhere.  In  particular,  the  radius  of  curvature  at  the  stagnation  point  is 


M  \W- 

•»  =  *.«’«  I 


[<’-2] 


Then  A„  at  the  stagnation  point  may  be  evaluated  by  using  relations  |-I4]  and  | ({:)].  We 
obtain 


Therefore  the  criterion  for  either  the  validity  of  the  IturUe-Schumann  iissuinption  or  for 
extinction  is  practically  independent  of  the  transport  properties.  They  enter  only  through 
the  values  of  Pr  and  Sc. 


